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B physics has played a prominent role in investigations of new physics effects at low-energies. 
Presently, the largest discrepancy between a standard model prediction and experimental measure- 
ments appears in the branching ratio of the charged current mediated B — > tVt decay, where the 
large r mass lifts the helicity suppression arising in leptonic B decays. Less significant systematic 
deviations are also observed in the semileptonic B — > D^*'TVr rates. Due to the rich spin structure 
of the final state, the decay mode B — >■ D'ti/t offers a number of tests of such possible standard 
model deviations. We investigate the most general set of lowest dimensional effective operators 
leading to helicity suppressed modifications of & — >■ c (semi) leptonic transitions. We explore such 
contributions to the B — >■ D*Ti>T decay amplitudes by determining the differential decay rate, lon- 
gitudinal D* polarization fraction, D* — r opening angle asymmetry and the r helicity asymmetry. 
We identify the size of possible new physics contributions to these observables constrained by the 
present B — ^ D^*'TVr rate measurements and find significant modifications are still possible in all 
i^ of them. In particular, the opening angle asymmetry can be shifted by almost 30%, relative to the 

standard model prediction, while the r helicity asymmetry can still deviate by as much as 80%. 
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Oh 
'^ I. INTRODUCTION 

(N 

^_^ Several recent experimental results in B physics have significantly constrained the possibility of large New Physics 

f~| (NP) effects in rare B processes. In particular, new CP violating effects in the Bs — > J I4'4> decay are already 

Q^ constrained to be of the order of the Standard Model (SM) expectations [1]. Similarly, the recent LHCb bound on 

JL. Br{Bs — > ^'^jjT) [2] already rules out any significant enhancement over the SM prediction in this decay. 

Qj On the other hand, existing measurements of the branching ratio for the charged current mediated B -^ tVt process 

r^ yield results which are systematically higher than the SM expectations [3] and the current world average for this lepton 

'~~' helicity suppressed decay rate is a factor of 2 above the SM predicted value [4] . B meson decays with r leptons in the 

— 1 final state offer possibilities of significant NP contributions not present in processes with light leptons. Namely, the 

^ large tau mass can uplift the helicity suppression of certain (scmi)leptonic decay amplitudes which are unobservable 

■^ in decays with light leptons in the final state. NP models with extended ElectroWeak (EW) symmetry breaking 

lO sectors - the Two Higgs Doublet Models (THDMs) being the canonical examples - often predict enhancements in 

^^ such helicity suppressed amplitudes. Existing studies of the B —> DtVt decay [5-8] have already shown how such NP 

^^ effects can be over constrained, and how additional kinematical observables in the three body decay offer tests of the 

C^ underlying short distance contributions not possible in the two body leptonic mode. 

<^^ In this respect, the B -^ D*TVr decay having two detectible particles of non-zero spin in the final state {D* , r) offers 

^N the opportunity of an even more complete investigation of the structure of possible NP contributions to 6 — > ctVt 

\ [ transitions [5]. The experimental reconstruction of the D* in the D-k final state allows to obtain the helicity structure 

>■ of this state directly. Similarly, the r lepton helicity can be inferred from its decays to tti^t. final states [6, 7]. This means 

that a number of experimental observables sensitive to possible NP effects, can be introduced. In the present study, 

, . we explore several such observables, like the differential distribution over the lepton invariant mass, the longitudinal 

j^ D* branching fraction, a D* — r opening angle asymmetry, as well as the r helicity fractions. To this end we employ 

a model independent effective field theory approach and identify NP contributions, which naturally predict helicity 

suppressed contributions in (semi) leptonic B meson decays while preserving the well established SM form of charged 

lepton currents (c.f. [9] for a recent related discussion). 

The paper is organized as follows: in Sec. II we investigate NP inducing helicity suppressed contributions to 6 — >■ c(.D(^ 
within the effective theory approach and evaluate existing constraints coming from the B — )■ DtD-t rate measurements. 
In Sec. Ill we focus on the various kinematical distributions and spin observables accessible in the B — >■ D*Ti>T mode 
and estimate their sensitivity to such NP, while the explicit derivation of the relevant polarized differential rates and 
evaluation of the corresponding hadronic matrix elements is relegated to the Appendices. Finally, we conclude in 
Sec. IV. 



II. NEW PHYSICS IN b ^ crv-r 

Following [9], we consider effective weak Hamiltonian, relevant for b — > clvi transitions in presence of NP con- 
tributing only to charged current interactions of quarks, while manifiestly preserving the well tested universal V ~ A 
structure of leptonic charged currents 

Heff =^^^ JbCM E (^V^L^^)+h.C., (1) 

where Pl,r. = (1 T 75)/2, while J^ is the effective 6 — )• c charged current. In particular, we are interested in 
NP contributions, which lead to charged lepton helicity suppression in B — )■ D^*^iPg, and are thus inaccessible in 
semileptonic decays to light leptons [9]. In general, this is the case if the NP contributions to J^ can be written as 
a total derivative of a scalar operator. In the effective field theory expansion, the most relevant (lowest dimensional) 
contributions to Jj^^ are then 

Jbc = c^PLh + gsLid^ (cPLb) + gsRid^ (cPRb) , (2) 

where the first term corresponds to the SM charged current, while gsL,SR sue dimensionful NP couplings. If the NP 
contributions are associated with a high NP scale Anp ^ fEW) then gsL,SR ~ 1/Anp- A particular and well known 
realization of such NP contributions is the THDM type II where only gsR receives a significant contribution. It is of 
the form gsR^ ^ —rrib tan^ P/rn^j^ [10] where tan/3 is the ratio of the two EW condensates in the model, and r7i^+ is 
the mass of the physical charged Higgs boson. Further NP contributions to J^^ relevant for helicity suppressed decays 
can be obtained via insertions of the d^ operator and are thus necessarily suppressed by at least two additional powers 

of Anp- 

In specific models one may relate NP effects inh -^ c transitions to other sectors, i.e. &, s — > w (B — >■ tDt, K — >• ^v^j) 
or c ^ d,s {Di^g-^ -^ dui), c.f. [11, 12], resulting in a more constrained parameter space. In the present study we 
will however not assume any underlying flavor structure and focus exclusively on observables in the 6 — >■ c sector. 
Before exploring such NP effects in various kinematical distributions of the B — >■ D*tDt, we need to consider existing 
constraints coming from the measurement of the B — >■ Dtv-t decay rate. In particular it turns out that most hadronic 
and SM parametric uncertainties cancel in the ratio between the tau and light lepton branching ratios [8, 13], i.e. 

^ Br{B ^ DrPr) , . 

~ Br{B -^ DeDe) ' ^ ' 

This ratio can already be predicted with considerable precision in the SM, and present estimates using either Lattice 
results with i^g^'' — 0.296(16) (updated value based on [8] using the recent precise form factor shape determination 
from [14]) or heavy quark expansion with -Rgjyf = 0.302(15) [7] agree within the stated errors. NP of the form in (2) 
results in a modification of the R ratio between the tau and light lepton rates 

R/RsM = 1 + 1.5Re[mr{gsR + gsL)] + l.O\mr{gsR + 5sl)|^ , (4) 

where we have again updated the expression in [8] using the form factor shape determination from [14], and the MS 
values at the ms scale have been used for the bottom and charm quark masses. Comparing these expressions with 
the experimentally determined values [16] 

Br{B+ -^ D^T+Vr)o^p = (0.77 ± 0.25)% , 
Br(S+^i)Vi^f)exp = (2.23 ±0.11)%, for£ = e,/i, (5) 

one can obtain constraints ^ in the complex plane of {gsR + gsh) as shown in the left plot in Fig. 1. 

For completeness, we provide SM predictions for the branching fractions by using the experimentally measured 
decay rates to light leptons, i.e. assuming no NP in those modes. Using inputs from [14, 16] we obtain 

Br{B+ -^ D°T+Vr)sM = (0.66 ± 0.05)% , Br{B^ -^ D-T+Vr)sM = (0.64 ± 0.05)% . (6) 



^ The related neutral B decay modes have also been measured [16] in both D and D* final states, however with less significance compared 
to the charged B modes resulting in less stringent bounds. 
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Figure 1: The 68% (in darker green shade) and 95% (in lighter yellow shade) C.L. allowed regions in the complex plane of 
NP parameters appearing in the effective 6 — >■ c charged current in Eq. (2). Shown are fixed combinations (the orthogonal 
combinations are set to zero): gsR + gsL (left plot) bounded by B — >■ DtVt, Qsr ~ 9Sl (center plot) constrained only by 
B — > D*TU-r, and gsL (right plot) contributing to both modes. The best fit NP benchmark point in the later case is marked 
with the symbol *. 

III. B -^ D*TVr DIFFERENTIAL DECAY RATES 

We consider the decay of a i? meson to a polarized D* (of helicity +, — or 0), a r lepton of a given helicity 
(At = ±1/2), and Vr (with helicity A^ = 1/2) as mediated by the effective Hamiltonian of the form in Eq. (1). The 
relevant kinematical variables describing the three-body decay are q^ = {pb — Pd*Y, where pb,d* are the B and D* 
momenta, respectively, and the angle 9 between the D* and r three-momenta in the r — v^ restframe. The detailed 
derivation of the polarized double differential rates is given in the Appendix A with the final result in Eq. (A14). 

We focus first on the decay distributions in absence of tau helicity information. Summing over both tan helicities 
At- in (A14), we obtain 
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where |p | is defined in Eq. (A2), and Hmn are the relevant (g^ dependent) helicity amplitudes, defined in Appendix B. 
Performing integration over dcosO in (7), we obtain 
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in agreement with the well known result [17, 18]. The presence of NP quark charged currents defined in (2) only 
affects the Hot helicity amplitude and can be encoded compactly as 
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(9) 



In the numerical evaluation of such NP effects we use the MS values for the bottom and charm quark masses at the rrib 
scale. The task of extracting information on NP from the differential decay rates thus reduces to obtaining sensitivity 
to the Hot helicity amplitude. 

We start with the B — > D*tDt branching fractions by integrating (8) over q^. These are also the only already 
measured observables in these modes. As in the case of the B — >■ DtDt decay, most theoretical uncertainties, related 
to the evaluation of the hadronic form factors defined in Appendix B are significantly reduced if one normalizes the 
B — > D*TDr rates to the modes with the light charged leptons in the final state [13] - one considers the ratio 
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Figure 2: The difFerential ratios R*{q^) (left plot) and R*l{(1^) (right plot) as functions of q^ . The black dashed curves are the 
SM predictions, while the blue dotted lines denote predictions for the NP benchmark point (see text for details). The 95% 
C.L. allowed regions (due to existing constraints from the B — >■ D^*'tVt branching ratio measurements) for fixed NP parameter 
combinations (the orthogonal combinations being set to zero) gsn — gsL and gsL are shown in lighter green and darker red 
shades, respectively. The gsR case is almost degenerate with gsL and is therefore not shown. 



In this way we obtain 



R* =K*s^{l + 0.12RcK(.9sfl - 5sl)] + 0.05|m,(gsfl - 9sl)?] 



(11) 



where we find for the SM prediction (using the recent precise experimental extraction [22] of the relevant form factor 
ratios)^ 



Rly, = 0.252(3) , 



(12) 



The stated hadronic uncertainty is dominated by the estimate of higher order perturbative and power corrections to 
the heavy quark limit of the Aq/Ai form factor ratio which presently cannot be extracted directly from data (see 
Appendix B). At this level of precision, EM corrections affecting B — > D*eDe and B — > D*TVr differently could become 
important [19] but the related uncertainty due to such effects depends on the particular experimental setup and would 
require a dedicated study beyond the scope of the present paper. The above expressions are to be compared with the 
experimentally determined branching fractions [16] 
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From this we can again obtain constraints in the complex plane of {gsL ~ 9sr) as shown in the central plot in Fig. 1. 
We observe that while certainly being complementary to the B — > DtDt mode, NP contributions to the integrated 
B — > D*TDr branching fraction are much more diluted. 

We also note that at present, the experimental measurements are systematically above SM predictions in B — > 
D^*^TDr decays. This is clearly demonstrated in the special case, where we set gsR = and study the combined 
constraints from both decay modes in the right plot in Fig. 1. We observe that such a fit mildly prefers a non-SM 
solution with gsL — — 0.9GeV~^. In the following we will use this benchmark point to evaluate the discriminating 
power of the various B — > D*ti't observables. 

Again for completeness, we provide SM predictions for the B -> D*tVt branching fractions by normalizing the 
value of i?sM to the experimentally measured decay rates to light leptons, i.e. assuming no NP in those modes. Using 
inputs from [16, 22] we obtain 



Br{B+ -^ D*°T+Vr)sM = (1.43 ± 0.05)% , 



Br{B° 



>r)sM = (1-29 ±0.06)% 



(14) 



in agreement with previous estimates [13]. 

Let us next explore the influence of NP on the ratio of singly differential decay rates with tau and light leptons in 
the final state 
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The value is obtained by averaging over the B^ and S" modes, for which in absence of EM corrections, Rgj^ differs by less than 0.001. 




Figure 3: The differential asymmetries Ag{q'^) (left plot) and A\{q^) (right plot) as functions oiq^. The black dashed curves are 
the SM predictions, while the blue dotted lines denote predictions for the NP benchmark point (see text for details). The 95% 
C.L. allowed regions (due to existing constraints from the B — >■ D^*'tVt branching ratio measurements) for fixed NP parameter 
combinations (the orthogonal combinations being set to zero) gsR — gsL and gsL are shown in lighter green and darker red 
shades, respectively. The gsR case is almost degenerate with gsL and is therefore not shown. 

Here, dV^/dq^ is the differential decay rate to a light charged lepton, where helicity suppressed effects are negligible. 
Taking into account present constraints on the NP contributing in linear combinations gsL a-nd gsR — gsL we obtain 
the 95% C.L. allowed bands in the left plot in Fig. 2. We observe that significant effects in R*{q^) (and consequently 
R*) are still possible, especially if NP contributions are aligned with the gsR — gsL direction - if they appear in the 
form of a pseudo-scalar density operator to which B — >■ DtVt has no sensitivity. 



A. Longitudinal D* Polarization and the Opening Angle Asymmetry 



Since NP of the form (2) only contributes to longitudinally polarized D* {D*^) in the final state, an increased 
sensitivity can be expected by using information on the polarization of the D* , which can be inferred from the 
angular distributions of its decay products (i.e. D-k). In (8) only Hqq and H^t contribute I?2'Si leading to a prediction 
for the longitudinal rate, again normalized to the light lepton mode 



Kl^ 



_ Br{B -^ DlTVr) 
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= 0.115(2) {1 + 0.27Re[mr{gsR - 5sl)] + 0.10|m,(gsfl - .9sl)|'} 



(16) 



where we have also given the estimated hadronic uncertainty of the SM prediction. In addition to this inclusive 
observable, one can also study the singly differential longitudinal rate ratio R*]^{q^) defined analogously to R*{q^) in 
Eq. (15). The presently allowed ranges for this observable are shown in the right plot in Fig. 2. Compared to R*{q^) 
this observable clearly exhibits an increased sensitivity to NP contributions. 

While the D* polarization information can be extracted directly from the angular distribution in (7), this requires 
experimental fits to two-dimensional decay distributions. In order to avoid such challenges, we propose a simple 
angular (opening angle) asymmetry defined as the difference between partial rates where the angle 9 between the D* 



and T three-momenta in the t 



Aeiq^) 



T rest-frame is greater or smaller than 7r/2 
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In the decay modes with light leptons, this asymmetry (Ag) can be used to probe for the presence of right-handed 
6 — > c currents, since these contribute with opposite sign to H±± relative to the SM. In the tau modes, it is sensitive 
only to the real part of NP ggi^ — gsR contributions and thus provides complementary information compared to the 
total rate (or R*). The presently allowed ranges for the Ag asymmetry are shown in the left plot in Fig 3. We observe 
that significant deviations from the SM prediction in this observable are still allowed. Also note that in the SM this 
observable exhibits a zero crossing at q^ ~ 5.6 GcV'^ , while this is not necessarily the case in presence of NP. On the 
other hand, the inclusive asymmetry Ag integrated over q^ is very small in the SM with AggM = —6.0(8)%; for our 
NP benchmark point we obtain A^^np = 3.4% , but even values as low as —30% are still allowed. 
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Figure 4: Regions allowed by potential future 10% relative precision measurement of R\ (shaded in gray and bounded by dot- 
dashed lines), 10% precision determination of Ax (shaded in green and bounded by dashed lines), and 5% precision measurement 
of Ae (shaded in red and bounded by dotted lines) in the complex plane of NP parameter gsL appearing in the effective b -^ c 
charged current in Eq. (2) {gsR is set to zero identically). All observables are assumed to be SM like. The present B — )■ D^ Wvr 
decay rates' best fit NP benchmark point is marked with the symbol *. 

B. Using r helicity 

It has been pointed out recently [7], that the spin of the tau lepton originating from semileptonic B decays can 
be inferred using the distinctive tau decay patterns. Therefore it is beneficial to consider the B — >■ D*tVt decays 
with taus in a given helicity state (A,- = ±1/2). In particular, assuming the standard V — A structure of the leptonic 
charged current entering the relevant effective weak Hamiltonian (1), the Ar = 1/2 helicity final states are suppressed 
by the tau lepton mass. Using the derivation of the polarized differential decay rates in Appendix A we obtain 
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which has the explicit form 
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The presently allowed ranges for this asymmetry are shown in the right plot in Fig 3. We observe that also in this 
observable significant deviations from SM predictions can be expected. Even in the inclusive asymmetry, integrated 
over q"^, where the SM predicts A^sm = 0.829(15), our NP benchmark point yields A^np = 0.36, while even slightly 
negative values are still possible. 



IV. CONCLUSIONS 



Within the effective field theory approach we have studied the most general lowest dimensional contributions to 
helicity suppressed (semi)leptonic 6 — >■ c transitions and found that a precise study of the exclusive decay mode 



B — >■ D*tDt could clarify the possible existence of such non-SM physics. 

Most hadronic inputs entering the theoretical predictions for the B — >■ D*tVt- rates can be substantially reduced 
by normalizing to the B — > D*ei?g^ mode. Using the helicity amplitude formalism we have investigated the presence 
of r mass suppressed helicity amplitudes not observable in S — > D*£De with i = e, ii. We have estimated these 
contributions using heavy quark expansion including leading perturbative and power corrections and derived precise 
predictions for the B — > D*tVt branching fractions in the SM. In order to further refine these estimates, Lattice QCD 
results for the form factor ratio Aq/Ai (defined in the Appendix B) would be most valuable. 

The B — >■ D*TDr mode has an unique sensitivity to the pseudoscalar density operator which does not contribute 
to the B — > Dri'r decay mode, while the opposite is true for the scalar density operator. Consequently, the precise 
experimental study of both B — >■ D^*'tDt decay modes can be extremely useful in constraining these kinds of beyond 
SM physics, especially, since present branching fraction measurements of all these helicity suppressed semileptonic 
modes are systematically above SM predictions. 

Due to the rich spin structure of the B -> D*Th'r final state one can introduce new observables such as the 
longitudinal polarization fraction of the D* [R*]^), the D* — t opening angle asymmetry {A0) and the r helicity 
asymmetry {A\). The discriminating power of these observables is demonstrated in Fig. 4 where we show the 
constraints that such possible future measurements would impose on the NP parameter space for a specific case of 
gsR, for which present B -^ D'^*\i>r rate measurements prefer values away from zero (SM). We see that 0(10%) 
precision measurements of K^^^Ag.Ax could critically probe such effects. We have also determined the differential 
ratios i?*, i?^ and asymmetries Aq, A\ as a function of the leptons' invariant mass squared. All observables still allow 
for significant modifications of the corresponding SM predictions and could thus help to disentangle the short distance 
contributions to B — ^ D^*'tPt rates. In term, our study points out the importance of precision measurements of 
B — >• D*tVt at the LHCb experiment, as well as at the future Super B factories. 
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Appendix A: Kinematics and Helicity Structure of the B -^ D*TU-r 

Following [17] we consider the kinematics of the decay of a -B meson to a polarized D* together with a charged 
lepton - antineutrino pair. As throughout the paper, we will assume that the decay is mediated by interactions of 
the form (1), i.e. that the leptons are produced via the standard left-handed charged current as in the SM, while NP 
could modify quark charged currents. However, generalization of our results to non-standard forms of the leptonic 
charged currents is straightforward. We denote the momenta of B,D*,i,i' with pbtPd* ^ki^^ky respectively, while 
q = Pb — Pd" ~ kg + k^. We also introduce the angles 9 between the D* and r three-momenta in the t — Dt rest 
frame, as well as Xi between the plane of the charged lepton and antineutrino momenta, and the decay plane of the 
D* . Helicity basis vectors of the D* (vector) meson are denoted as e", while assuming standard lepton interactions, 
we can use e^ for the four basis vectors describing the total helicity of the charged lepton-neutrino system. In the 
rest frame of the B meson with z axis along the trajectory of the D* , a suitable basis for the lepton pair helicities is 



£^(±) = ^(0,±l,-z,0) 



1 

£p(0) = ^(|p|,0,0,-go), 

Vr 

e^(t) = ^(9o,0,0,-|p|), (Al) 



where go = {itt-'b ^ "*Z5* + '?^)/2mB and 



Ai/2(m|,m|).,g2) 
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with A(a, b, c) — a^ + b'^ + c^ — 2{ab + bc + ca). They satisfy the following normalization and completeness relations 

e*(TO)e''(TO') =5mm', for (m,TO' = t,±, 0) , (A3) 



^ £^(m)e*(m')5™„/ = 5^^ . (A4) 

ra,m' 

Similarly a convenient helicity basis for the D* is 

£„(±) = T-^(0,l,±i,0), 

eo(0)-— (|p|,0,0,-Bd*), (A5) 

where Eo* = {m\+m?jj, —q^)/2mB is the energy of I?* in the B rest frame. These basis vectors satisfy the following 
normalization 

£*(m)£"(m') = -^™™', (A6) 

and completeness relation 

El \ ( i\!r , PD'aPD'13 /.„N 



We can now introduce helicity amplitudes, H±±, Hqq and Hot describing the decay of a pseudo-scalar meson into the 
three helicity states of a vector meson and four helicity states of the leptonic pair 

Hmm{q^) = £{mY*H^,{m), for m = 0, ± , 
Hot{q^)^e{m^tY*H^{n = Q). (A8) 

Here, Hf^{m) is a corresponding hadronic matrix element, and m,n denote helicity projections of the D* meson and 
the leptonic pair in the B rest frame. 

If mediated by the effective Hamiltonian of the form (1) for arbitrary quark charged current J^^, the B — >■ D*tv^ 
triply differential decay rate can be written as 



(fT, ^ GllKfcp IpI / mj\ 
dq^dcosOdx (27r)4 2ml \ ^V ^"^ 
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where Lf^,j^Hf^i, are the leptonic and hadronic current tensors. Using completeness relations of the helicity basis vectors 
we can rewrite L^,yH^^ as 



mm' ,niv 



Following [20, 21] we can expand the leptonic tensor in terms of a complete set of Wigner's d'^ functions, reducing 
L^uH^'^ to the following compact form 

Xe,XD,,\t„,yf„:JJ' (All) 

where J and J' run over 1 and 0. In term, the lepton helicity amplitudes, hrxi,\^) for ^ left-handed weak current are 
given by 

kxi-x.) = ^M^ih'^i'^ - 7')«.(A.)?^(A,,), (A12) 

where for massless right-handed antineutrinos Xi, — 1/2 and Xi^ — Xe ~ Xi, in the £v center of mass frame by angular 
momentum conservation. It follows that the two non-vanishing |/i(Af,A„) P contributions are 

2 
\h-i/2,i/2\' = 2iq'-mj) and \h,;2,i/2\^ ^ 2^{q^ ~ mj) . (A13) 
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Finally, using the standard [16] convention for Wigner's d-functions and performing the trivial integration over x ^^ 
obtain 

SL (^^ = -1/2) = ^'^'S'"' f 1 - # V [(1 - cos ^)'^++ + (1 + cos (^fHl- + 2 sin^ ^^oo] , 



dq^d cos 9 25677^™^ \ q^ 



S^(^' = '/^) = %Sf (' - f)' F ['^'°'''^- + "--' + ^"'»' - '^"""^^'^l ■ '*"' 



dq'^d cos 6 25671^771^ \ q^ J q 

from which Eqs. (7), (8) and (18) can be easily derived via summation over A^ and/or integration over cos0. 

Appendix B: Helicity amplitudes and hadronic matrix elements 

In the SM (in presence of only the first term in (2)), the helicity amplitudes Hmn can be writen as 
Hftiq') ^ {ms+mn^)A,{q')T ^""^ |p|^(g^), 



Hl^iq' 



2m£)* y'q'' 



(m| - ml, q^){mB + mo^^iq^) - ^"^^'P'' A^{q^) 

niB + ruD' 



H-^{q^)=^-:^^^A,{q^). (Bl) 

vr 

where the form factors parametrizing the relevant hadronic matrix elements are defined as 

{D*{pD',e^)\c^^h\B{pB)) - "^'^y^ £^.affe*>gpg,. , (B2a) 

niB + rao* 

{D*{pD',ea)\c-i^,^5b\B{pB)) = 2mD' Aoiq'^) — ^g^ + (tob + "^z^*) ^i(9^) ( e^ ^Qt^ 

-A,{q') '*_■'' ({pB +PD^), "^T^^ g,) . (B2b) 

mB + mn* \ Q J 

In presence of NP of the form (2), one needs to evaluate two additional matrix elements given by 

{D*{pD',ea,)\cb\B{pB)) - 0, (B3a) 

{D*{pn-.,eo.)\n,b\B{pB)) = ] q''{D*\-cj^j,b\B°) = ^"^f' Ao{q^)e* ■ q . (B3b) 

TOfc + Trie ITT-b + ITT-c 

The final effect of such contributions to the differential rates can be encoded into the Hot helicity amplitude as given 
in Eq. (9). 

We can use further information on the form factors given by precise differential decay rate measurements in _B — > 
D*ivt [22], as well as their perturbatively computable properties and relations in the heavy quark limit for the &, c 
quarks [23, 24]. In this limit it is customary to employ a new kinematical variable 

W = VB-VD'^—^ — —. (B4) 

with w^, and u^, being the four- velocities of the B and D* meson respectively. One can then define an universal 
form factor 

HaAw) = A,{q^)-^ ^, (B5) 



and ratios i?i , R2 and i?o in terms of which 



2. _ Ro{w). 



A,{q')^^^hAAw) 



Rb* 

2, _ R2{w) 



A,iq') = ^^hAAw) 



Rb> 

2s _ Rl{w) 



V{q') = ^^hAAw), (B6) 
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where Rd' — '^^/mBrriD* jirriB + mo*)- The w dependence of these quantities in the heavy quark hmit reads 

HaAw) = /iAi(l)[l-8p2z + (53/92-15)z2- (231p2-91)z3] , 

R^{w) = i?i(l)-0.12(w-l)+0.05(u;- l)^ 

R2{w) = i?2(l)+0.11(u;-l) -0.06(u;- l)^ 

Roiw) = i?o(l)-0.11(w-l)+0.01(u;-l)^ (B7) 

where z = {y/w + 1 ~ ^/2)/{\/w + 1 + ^/2). The first three expressions can be found in [25], while we have derived 
the fourth using the resuhs of [25]. Above relations contain free parameters /i^j(l), i?i(l), i?2(l),p^, which can be 
extracted from the well measured B -^ D*£i>i decay distributions. In our numerical evaluation of the B -^ D*tVt 
differential decay rates we employ the results of a recent Belle analysis [22] . The virtue of this approach is that most of 
the associated hadronic uncertainties actually cancel in ratios of decay rates to tau versus light Icptons, as previously 
demonstrated for the case of B — > Drvr [8] . 

In addition to these inputs, the B -^ D*tDt rate also depends on Ro{l), which cannot be extracted from B — >■ D*£Dg 
studies, since it only appears in the helicity suppressed amplitude Hot- In the exact heavy quark limit [i?o(l)]HQET — 1- 
Leading order perturbative (in as) and power (l/rrii,,,) corrections are known [23, 24] for the linear combination 

i?3(l)^«fc^:l±*g^^i±^:l^=0.97. (B8) 

(1 — r)^ 

The same calculation predicts i?2(l) — 0.80, yielding eventually i?o(l) — 1-22 . Experimentally [22] however, i?2(l) = 
0.864(25) and inserting this value into Eq. (B8) yields our final result -Ro(l) — 1.14. In all our numerical calculations we 
use directly Eq. (B8) and conservatively assign a 10% uncertainty to this value accounting for higher order corrections. 
In the future, a more reliable and precise determination of -Ro(l) (or equivalently i?3(l)) could be obtained on the 
lattice, similarly as has already been done for the helicity suppressed B ^ D matrix elements [15]. 
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